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Abstract 

Electron transport phenomena in disordered electron systems with spin-orbit coupling in two dimensions and below are 
studied numerically. The scaling hypothesis is checked by analyzing the scaling of the quasi-lD localization length. A 
logarithmic increase of the mean conductance is also confirmed. These support the theoretical prediction that the two 
dimensional metal in systems with spin-orbit coupling has a perfect conductivity. Transport through a Sierpinski carpet is 
also reported. 
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1. Introduction 

Research in spintronics and carbon nanotubes 
sparked renewed interests in the effects of spin-orbit 
coupling on quantum transport [1,2]. One of the re- 
markable features of non-interacting disordered elec- 
tron systems with spin-orbit coupling is the existence 
of a metallic phase in two dimensions (2D) [3,4]. More 
exactly the 2D metallic phase appears in systems with 
symplectic symmetry, i.e., in systems with time re- 
versal symmetry but without spin rotation symmetry. 
Such systems are an exception to the prediction of 
Abrahams et al. that there is no metallic phase in 2D 
disordered electron systems [5] . 
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Recently some 2D metals have been discovered ex- 
perimentally [6] . Although the mechanism inducing the 
2D metal might be different from spin-orbit coupling, 
we think that detailed understanding of the 2D sym- 
plectic systems gives some insights on general trans- 
port properties of 2D metals. 

A scaling theory for the symplectic systems [3,7,8] 
can be described by the (3 function for the conductance 
g (in units of e 2 /h) 



ding 
dlnL' 



(1) 



where L is the linear size of a system. (The conductance 
is actually a distributed quantity. See Refs. [9,10,11] 
for discussions on the scaling hypothesis of the conduc- 
tance.) Because of the so-called anti-localization effect, 
the P function shows non-monotonic behavior in the 
symplectic systems. The asymptotic behavior in the 
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strongly metallic region in 2D is conjectured to be 

/3(ln 9 ) = — (2) 

This indicates that there is a metallic phase in 2D. 

In the scaling theory for the symplectic systems in 
2D, however, there is a long standing problem for more 
than two decades [12,13]: the logarithmic divergence of 
the conductance in the 2D metallic phase in the limit 
L — ► oo indicated by Eq. (2). Since the conductivity 
and the conductance are the same quantity in 2D, this 
predicts a perfect conductivity, a = oo, in spite of the 
system being disordered. This is in contrast to 3D dis- 
ordered metals, in which the conductivity converges to 
a finite value. 

According to the renormalization group theory of 
continuous phase transitions, the scaling hypothesis is 
valid when the correlation length is much larger than all 
other microscopic lengths. The scaling argument which 
deduced (2), however, is based on the weak localization 
theory in which disorder is supposed to be weak. In such 
a weakly disordered metallic phase this condition on 
the correlation length may not be satisfied. Therefore, 
a numerical check of the scaling hypothesis in the 2D 
metallic region is desirable. 

In this paper, we report a numerical check of the 
scaling hypothesis in the 2D metallic region. We also 
report numerical calculation of the Landauer conduc- 
tance in 2D, which is a direct check of the prediction 
of a perfect conductivity. Last, numerical simulations 
of transport through a Sierpinski carpet are reported. 



2. The SU(2) model 

Among the models with symplectic symmetry, we 
employ the SU(2) model [14,15], 

H = £i4 a Civ - R{i,3) aa ,c\ a c i<yl . (3) 

The random potential e» is distributed with box distri- 
bution in the range [—W/2,W/2]. Random spin-orbit 
coupling is included in the nearest neighbor hoping 
term. We parameterize the hopping matrix as 



R(i,3) = 



e' a ^ cos Pi- 



sin 



-e i7lJ ' sin fa e cos/3 i: 



(4) 



and we distribute ctij and jij with uniform probability 
in the range [0, 27r) , and (3ij according to the probability 
density p(p)df3 = sin(2/3)d/3 in the range [0,7r/2]. In 
the actual simulations R(i,j)'s in one direction are set 
to the unit matrix with the aid of the local SU(2) gauge 
transformation. 

The SU(2) model has an advantage that corrections 
to scaling arising from irrelevant variables are smaller 
than other models [14,15], so it is quite a useful model 
when we study universal aspects of Anderson localiza- 
tion and the Anderson transition. 



3. Numerical check of the scaling hypothesis in 
the 2D metallic phase 

We analyze the quasi- ID localization length A on a 
quasi-lD strip with width L [16,17]. Periodic bound- 
ary conditions are imposed in the transverse direction. 
The scaling hypothesis implies that the dimensionless 
quantity A = X/L obeys 



A 



(5) 



where £ is the 2D correlation (localization) length 
which depends on W and the Fermi energy E. The sub- 
script ± indicates the metallic and insulating phases. 

We have calculated A for sizes L — [16, 128] with 
an accuracy from 0.3% to 1.0%. The ensemble transfer 
matrix method has been used [18]. In addition to the 
data at E = 1 analyzed in Ref. [15], we have also ac- 
cumulated data at E — 2 in W = [0.0, 4.5]. From the 
numerical data the correlation length £ at each pair of 
(E, W) and the scaling function F+(L/£) are estimated 
using numerical fit described in Ref. [15]. The quality 
of the fit is assessed by the goodness of fit probability 
Q. To eliminate the ambiguity of the absolute scale of 
£, we set £ = 10 at (E, W) = (1, 3). 

Figure 1 demonstrates the single parameter scaling 
in the metallic region. All data fall on a common scaling 
curve when A is plotted as a function of L/£, indicating 
the validity of the scaling hypothesis. In a strongly 
metallic region (A ^> 1), the data are well fitted to 



A = 6 + cln(L/£) (A>1), 



(6) 



as shown in Fig. 2. 
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Fig. 1. Single parameter scaling in the metallic region in 
2D. The goodness of fit Q is 0.3. 
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Fig. 2. Single parameter scaling in the strongly metallic 
region in 2D. The goodness of fit Q is 0.7. 



In Ref. [15], the (3 function for A from the metallic to 
the insulating limits was estimated. The estimate has 
not changed significantly when we add data for E — 2. 



4. Logarithmic increase of the conductance in 
the 2D metallic phase 

In last section and in Ref. [15] , the scaling hypothesis 
has been checked numerically. This supports the scal- 
ing theory for the 2D symplectic systems. Furthermore, 
the logarithmic increase of A in the strongly metallic 
region gives the impression that the conductance may 



Fig. 3. The mean conductance (g) vs. system size L in the 
2D metallic region. The solid lines are a guide to the eye 
only. For a reference, (g) = a + tt~ 1 InL is indicated. 



also increase logarithmically. In this section, we study 
the size dependence of the mean conductance directly. 

We now attach two perfect leads to the square sam- 
ple with linear size L and calculate the two terminal 
Landauer conductance g (in units of e 2 /h) with the 
transfer matrix method [19]. A fixed boundary condi- 
tion is imposed in the direction transverse to the cur- 
rent. The hopping matrices in the transverse direction 
are set to the unit matrix, and the Fermi energy E = 1. 
We have accumulated 10000 samples for L = [64, 256] 
and 5000 samples for L = 384. 

As shown in Fig. 3, the mean conductance (g) in 
the strongly metallic region does show the logarithmic 
increase 

(g) & a + 7T -1 InL «ff)>l), (7) 

in agreement with the prediction of the scaling the- 
ory [3]. The prefactor of the logarithmic term, which is 
an important parameter for the scaling in the strongly 
metallic region, is also consistent with it. Note that a 
perfect conductivity does not mean perfect transmis- 
sion. The transmission probability per channel rather 
goes to zero as L increases, so most of incoming elec- 
trons are reflected. 

It is left for future to understand physics behind this 
possible perfect conductivity in 2D. 



3 



4.0 
3.8 
3.6 
3.4 
3.2 
3.0 
2.8 
2.6 



y'a+7T ] \nL 




Vy - 


-U 


y' ' u- 


--2 



10 20 



50 

L 



100 200 



500 



Fig. 4. The mean conductance (g) vs. system size L = 5 k for 
the Sierpinski carpet SC(5, 1, k). The solid lines are a guide 
to the eye only. The dashed line indicates the logarithmic 
increase in 2D. 



5. Transport through a Sierpinski carpet 

Recently we have obtained numerical results which 
suggest that an Anderson transition occurs even be- 
low 2D in the presence of spin-orbit coupling [20] . This 
is based on numerical simulations of electrons on a 
Sierpinski carpet SC(5, 1, k), where k is the generation 
number. In the limit k — * oo, SC(5, 1, k) becomes a true 
fractal whose spectral dimension is d s = 1.940 ± 0.009 
[21]. 

Here we study the size dependence of the Landauer 
conductance through the Sierpinski carpet in a delo- 
calized phase. We have attached two perfect leads to 
SC(5, 1, k) and have calculated the conductance with 
the recursive Green's function method [22]. A fixed 
boundary condition is imposed in the transverse direc- 
tion. The hopping matrices in the transverse direction 
are set to the unit matrix. We set the Fermi energy 
E — 0.9 and the disorder W = 0, 2 where relatively 
large mean conductance is obtained. The number of 
samples is 10000 for k = 1, 2, 3 and 500 for k = 4. 

Figure 4 shows the mean conductance (g) as a func- 
tion of the linear size L = 5 fe . The mean conduc- 
tance increases with L, possibly indicating a delocal- 
ized phase. It also indicates that the increase is slower 
than that in 2D. It is an open problem whether the 
conductance in the possible delocalized phase of the 
Sierpinski carpet diverges. 
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